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ABSTRACT 


An attempt has been made to find out the Dnections Of Arrivals (DOA) of narrow band signals 
(which aie non coherent) leceived ovei a umfoim Imeai airay The DOA algorithms like Multiple 
Signal Glassification (MUSIC) are critically dependent on computation of noise eigen vectois The 
MUSIC algoiitlmi involves any of the matiix iterative procedures of the complex Heimitian matrix 
(auto coiidation matiix of the piocess) is bound to have the limitations posed by the paitiiulai 
maliix iteiative proceduie both in teims of computational time and complexity For fast leal 
tune applications alternatively Direct Noise Snbspace Basis (DNSB) method has been studied for 
1C il data DNSB method is simple rncl computationally efficient Unlike the above mothocls the 
Subsp ue Based Approach which diiectly uses infoimation from both the noise and signal snbspaces 
has been veiilied ( ompaiisions with the above tliiee methods aie made For one source case the 
( R botrnd his been deiived and veiihed 
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Chapter 1 

Introduction 


The piaclical pioblems of interest m aiiay signal processing aie extiactmg the desiied paiameters 
su(h as tli( diuctions of anival (DOA) powei levels and cioss coirehtions of the signals present 
in Ihc scene flora the avulable infoimation including the measured data Often one may also be 
specific ally inteiested m the actual signal of one of these sources and m this case it is necessary to 
estimate the actual waveform associated with the desired signal by improving the overall reception 
111 an cnviionraent having seveial sources To achieve this ideally it should be possible to suppress 
the undesiied signals and enhance Ihe desired signal 

Sensoi airays have been m use for several decades in many practical signal processing applica 
tions Snell an an ay consists of a set of sensors that aie spatially distiibuted at known locations 
with reference to a common refeience point These sensors collect signals from sources m then 
held of view Depending on the sensoi chaiacteiistics and the path of propagation the source 
wavefronts undergo clitermimstic and/oi random modifications The sensor outputs are composed 
of these souice components and adclitive noise such as measuiement and theimal noise [ 1 ] 

Ihe most common means for distinguishing a desired signal from an mteifexence signal is the 
existence of a known ficcpiency band within which the desiied signal must fall We have to know the 
desned signal well enough so that it can be distinguished from interference signals it is reasonably 
to expect that the nature of the desired signal is known even though certain signal paiameters 
(such as DOA, amplitude and phase) may have to be estimated 
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1 1 Problem Formulation 

Signal IS modelled as a sample function of a random process Assumptions are made that the 
transmission medium is homogeneous so that the signal wavefront is perfetly coherent over the 
entire array(i e the signal components received by various sensors are delayed rephcas of each 
other) The noise is considered additive white gaussian (AWG) and incoherent from sensor to 
sensor 

In the fig 1 1 the plane wavefront is making an angle 0 on the array broad side with the line 
joimng sensors in the hnear array, d is the distance between the sensors and it is normadised with 
respect to the half wavelength associated with narrow band signal y(t) Also dj are 

similarly normalised distances of the sensors from the reference point 

The problem of Directions of arrival(DOA) estimation is solved and analyzed using various al 
gonthms hke MUSIC(MUltiple Signal Clasification) method, Direct Noise Subspace Basis (DNSB) 
method, and Subspace based approach to parameter estimation(here. Directions of arrival) 

1 2 Organization Of Thesis 

Chapter 2 contains a concise but complete in itself description of the important Directions Of 
Arrival (DOA) estimation algorithms 
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( liiptei 3 deals with the case of single source m detail It gives the theriotical lower bound 
on vaiiance of the bais of the parameter a; = tt cos 0 which when infinite snapshots are taken then 
the ( R bound depends on the SNR and number of sensors Heie we assume that our estimator is 
unbiased 

Chapter 4 gives the simulation results followed by a breif discussion 



Chapter 2 

Dirctions Of Arrival Estimation 


Tho complex envelop ol the total received signal at sensor is given by [2] 

= y{t)eip{{-jwd cosd)) + n (t) (2 1) 

n,(i) IS the noise at the sensor and 

E[n^(t)n^(t)] = cr^S j (2 2) 

The genenlisation ol signal model at the i*''’ sensor foi M identical sensors receiving signals fiom 
A nniiow band signals 1/1 (i) j/2(0) i/jv ( t) that arrive at the array from directions O 2 Gr 
with icspect to line of array is given by 

A 

~ yr{t)e%p(—j{'n'd cos ^ (0 '^) 

=1 

Wheic aie the input sinusoidal fiequencies Dropping the sinusoidal t dependence m the 
summation we have 


A 

a,(t) = 'Y^yr{t)exp{-j-Kd cos dr) + 

7=1 

The aiiay output vectoi Y(l), therefore can be wiitten as 

X(t) = AY(t) + nit) 


(2 4) 


(2 5) 
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wheie i (L) = [yi{t) yi^t) 2/7v(0]* ( signal vector) 
n{t) = [ni{t) nzit) nM(.t)f (noise vectoi) 

I = \/M[au,j ] IS called the diiection vector matrix of size M X R where wa, = tt cos 0a 

and 

fl(wA) = diieclion vectoi (of size M) is given by 

exp(-jdMUJk)f (2 6 ) 

IS used foi noxmalisation of the diiection vector Under the assumption that the signals and 
noises ire stationaiy zoiomean iincori elated landoin piocesses the array output covariance matrix 
IS given by 


S = JC[X{t)\^{t)] = ARA^ + cr^I (27) 

R IS the source covanantc matrix ol size A x A 

In the iollowmg section we will briefly describe some of the important algoiithms that exist for 
the solution of the Diiection of aiiival(DOA) pioblem We can classify them as Eigen stiucure 
based methods and non eigen structuie based methods MUSIC PISARENKO ESPIRIT GEESE 
etc comes undei foimei one where as Bartlett s Beamformer method Capon s Minrmum variance 
estimator and Linear prediction method comes nndei latter case First we will desciibe the latter 
case 

2 1 Non Eigen Structure Based Methods 

2 11 Beainformer 

Bcaxiiloimmg ordinarily involves foiming multiple beams fiom multielement ariays through the use 
of ippropriate delay and weighting matrices The conventional beamformer introduces apprpiiale 
delays at the output of the each sensor to compensate the piopagation delays of the wavefiont 
reaching the array [3] As the name implies, the array output weights are chosen to be phase 
factors required to stem the an ay along some specific direction 9 i e 
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Figure 2 1 Beamforming 



Where R = E[X{t)X^{t)] 
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In a single taiget seen 1liis estimatoi measures the actual power when steered along the true 
diiection of aiiival lesultmg a single peak m that direction But it can be easily verified that if two 
souices aie located inside the mam beam the peaks m fact meige m to a single peak resulting in 
a loss of resolution 

2 12 Capon’s Minimum Variance Estimator 

In gcnenl the iiiay output powei contains power fiom the desiied signal along the look direction 
as well as power in undesned ones along othei direction of ariival To minimize the contributions of 
this Intel set the aiiny output power is minimized while keeping the gam along the look diieclion 
constant [4] 


min RW subject to |u/^a(a;)| = 1 


(2 12 ) 


The solution to the weight vectoi foi the positive covariance niatiices is 


substituting m ecpi (2 11) foi fl(m) 




R ^a(a;) 
a^^{uj)R~'^a{u)) 


(2 13) 


Pc(^) 


1 

a^(u})R~^a{uj) 


(2 14) 


The estimator m eepr (2 14) has supeiioi resolving power compared to the Bartlett beainformei 
output 

2 13 Linear Prediction Method 

In this case one of the sensoi output is predicted as a linear combination of the remaining (M - 1) 
sensoi outputs at my instant and the piedictor coefficients are selected so as to minimize the mean 
squire error 


Letting a, I i„_j 
for Xn we have 


Xn-M+i stand for the predictor for the M sensor outputs and Xn the predictor 
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M~\ 

r, =-^aXn_i (2 1^)) 

= 1 

Tins selves ilio eiior as 

M-\ 

fn “ ^ ^ d X.fi^ ao = l (216) 

=0 

Minimization of the mean squaie eirorjE|f with respect to the unknowns results in the following, 
standaicl set of lineai equations 

M-l 

-C[c, T a] = a jC[a n-» 2, = 0 A, = 1 2 M - J (2 17) 

1=0 

and the RnU mean sq^lare eiioi is given by S\e, ^ = 8m-i Undei the assumption that the rioss 
toirclation between the sensoi outputs (the sensoi outputs aie spatially wide sense stationaiity) 
depend only upon intei element distances then 

B[x, _ = r(^ - 0 = 

and the cqn (2 1^) becomes 

M-l 

^ a, r(A,-i) = 0 1.= 1,2, M -I 
1=0 

md also 

M-l 

^ a r{-i) = 8m-\ (2 20) 

t-O 

Thus Xn can be thought of as the output of a linear shift invariant system diivtn by an 
imconelated noise piocess with averege power 5iw-i as shown below 

If the enois e,i e, aie also nncorrelated for all n k then the input repiesents a wlute noise 
process and Xn an autoregressive process of older M — 1 Prom hg 2 2 the power spectial density 


(2 IS) 


(2 19) 
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Figure 2 2 Linear prediction model 
5 (cj) = through the relation 

S.{u) = \G(^^p(luf S,(u.) = (2 21) 

and 

H{z) = H{exp{ju>)) = 1 + ai 2 "^ + + (2 22) 

The above polynomial has its zeros close to the unit circle and gives sharp peaks m the output 
spectrum correspond to the DOA estimate of the signals present The mam disadvantage of the 
hnear predictor is that it exhibits spurious peaks in the output spectrum whenever the array 
pointing direction is away from one of the actual direction of arrival [2] 

2 2 Eigen Structure Based Methods 
2 2 1 MUSIC 

In the fig 2 3 shown below, there are M identical sensors and receiving signals from K narrow band 
signals yi{t) i/ 2 (f), VKii) that arrive at the array from directions 0 i, 02 , wrt the line of 

array 

We assume here that none of the signals are coherent i e , |p,j | ^ 1 Using superposition of 
signal contributions of all the incident planewaves, the received signal x,{t) at the t*'* sensor can be 


written as 
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Figure 2 3 A typical array scene 


k 

= E + n,it) (2 23) 

tel 

where is the sensor element position, normalized distance to half of the common wave 
length, wrt some common reference point, and n^(t) stands for the additive noise at the sensor 
It IS assumed here that signals and noises are stationary zero mean, uncorrelated random processes, 
and further 


E[n,(t)nf (t)] = 


From the eqn (2 5) the array output vector X{t) — i4y(t) + n(t) 
where Y {t) = [t/i(t), y 2 (t), , yh (t)]* 

n(t) = [ni(t),n 2 (<), ,njw(<)]* and 


1 

g-jwi 


1 


(2 24) 


(2 25) 
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IS sn (M X A ) m-xtnx consisting of A diiection vectors From our assumptions M x M array 
oiilput rovaii'incc matiix S' = C[X{l)\^{t]] Ins the form 

5 = AE[Y{t)Y^{i)] + E[n{t)n^it)] = ARA^ + a^I (2 26) 

Here we set that covariance matiix is lieiinitian and positive definite and souice covariance 
nntiix E = E\Y{t)Y^\t)\ of size A x A which is non singular when there are no coherent sources 
aie present in the data 

Assuming M > A and that A is full lanh A which implies that nonnegative definite matrix 
4 A A" is also of rank A Therefore S will have M - h lepeated lowest eigenvalues {=cr^) 

llieH foK if V is in eigenvectoi of S foi the eigenvalue cr^ we have 

ARA^^V = 0 (2 27) 

Unitary dngonalizability of S implies oithonormal decomposition of the data space in to noise 
and signal subspaces of dimension M - K and A respectively 
Eqn (2 27) can be further simplified to 


A^^V = Q (2 2R) 

St lied in words the eigenvectors associated with the repeating lowest eigenvalue of 5 are or 
Ihogonil to the direction vectors coiiesponding to the actual arrival of angle 

However owing to the presence of noise in the eigenvector estimates that span these two sain 
pie subspaces(signal and noise subspaces) the oithogonahty relations of eqn (2 28) are no longer 
valid The best that we can do is to search for the signal vectors that are most closely orthogonal 
to the noise subspace Accordingly in the MUSIC algorithm it is proposed to estimate the angu 
Hr fiequevcies of the complex sinusoids m the input signal as the peaks of the following sample 
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spocti um 


S{u) 


1 

E^=K+^\aHiu)V\^ 


(2 29) 


Heie we bcc that loi the above equation to be meanmgM M > K + 1 and this restricts the 
minimum numbei of tequiied sensoi elements to atlea&t one more than the total number of sources 
picseiit in I ho scene 


In this clnptei wo biieiiy explain Diiect Noise Subspace Basis method and Sub Space based 
appjo acli to est imale the diiections of aiiival from 5 the correlation matrix which we shall assume 
lus been finly well calculated (by possibly fan sample mean of data snap shot covariances) 

2 3 Direct Noise Subspace Basis Vectors Finding 

ARA^V = 0 if and only if A^^V = 0 


Now we shill develop an algoiithm to compute V directly from the correlation matiix which 
deteimmcs noise sub space Since the zeios of (2 28) give the DOA estimates it is necessary to 
find all the M - K vectors V satisfying (2 28) 


To simplify cleiivxtions we shall assume a uniform end clement lefeienced and seperation nor 
malised an ly so that 1, i = 1 2, , M Thus (A„ , ) = m = 1 M n = I K 

iml tlicic foie, {A,n „) = m = 1 A n = 1 M Also let P be the right shift operator 

oil jfM 

SO th it /' IS i matnx is 


0 0 O' 

1 0 0 

0 1 0 


(2 30) 


and = 0 


0 


1 


0 
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2 3 1 Theorem 1 

let, M > 2A And = [l|Ao||Bo] with Aq = A x A and Bq = (A x (M — A — 1)) be a partition 
of A Then Ihe set v Pv where v = [1|(Aq^ 1)|0 0]^ forms a set of basis for noise 

subsp ICC of 4 

For the pioof of this theoiem please lefer [5] 

C iveii Ihe above general form of the siibspace basis vectors we next attempt to calculate these 
fiom the data Since Ihese are the basis vectors of of 5*0 = 5 — cr^I we get from (2 26) that, 


Sol/ = ( S' - <jH)V = ARA^^V = 0 


(2 11 ) 


f c( us p u I ition So as follows 


S’! So 

S2 54 



0 


wheio S*! IS (A + 1) X (A + 1) and other S', ’s are compatible This gives 


51 

52 


rank 


1 

Vi 

Vk 


= 0 


51 

52 




Again Icl s be the hist column of 


51 

52 


(2 32) 


(2 33) 


(2 34) 


and use the partition 
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to give 



(2 15 ) 


(2 36) 


Ileio it cxu bo seen that laiil S'" < A (the eqiiality holds always if M < 2A — 1) In Ihe 
present CISC if S“ is full lank v is diiectly detoinimed once eqn (2 16) is solved by using the left 
invciso ol S~ 


Ihe 1 ink of 5“ may not satisfy the above eqn(2 34) in case of the data in real environment 
Then by using methods like Gnusian elimination or similar method of rank testing we can obtain 
good lesnlts 

2 4 Subspace Based Approach To Parameter Estimation 

llus xppioacli [b] &1 aits with a model of second order data(C orrelation matrix) geneiated fioiu tlie 
samples ofM tcceived signals '!t(n) ? = 1 2 M Each of these signals contains either a known 
01 eslimUed uuml)ei K of model xtely correlated signals Vkin) k = 1 2 R having distinct 
dilution ol u iivals(DOAS) di Qi Ok with lespect to a linear aiiay ol point sensors embedded 
m additive independent and idenlically distiibiited(i i d) Gaussian noise 

Ihc dimension of the signal subspace is estimated by comparing the magnitudes of the eigen 
values of the coii elation matrix and using a thieshold 

A 

fc=i 
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The ibove model provided in (2 37)with respect to the direction of ariival(DOA) pioblem appilies 
also in s( V( lal olhei situations including the problem of parameter estimation of sinusoids 

2 4 1 Analysis of the standard model 

The stand iid model foi the (M x M) con elation matrix 5 was given in (2 26) 

S = E[x{n)'T^^ (n)] = ARA^ + (2 18) 

Suite S IS Ileimiti in snd is modeled as in (2 38) it must be diagonalizable by a unitary matiix 
B Thciofott 


= B^({ARA^^) + cr^I)B 

= Dtag[(iii + a^) (/tA+o'^)cT^ cr^] (2 39) 

is a (M X M) diagonal matiix Foi notational bievity we deifne below the (A X A) and {{M — 
A) X [M — A)) diigonal matiices Ai and A 2 , lespectively 

A] = Diag[\iX'i A/^j A = cr^ for 2 = i A (2 40) 

A 2 ^ Dtag[XK^i Xm] \h-i- = for i=l M-K (-2 41) 

DtnoU the {M x A) matiix composed of the eigen vectois associated with the eigenvalues 
\i \2 \h by Bs (signal subspace oigenvectoi matiix) and 1he(A/ X (M — K)) matrix composed 
ol the (M - A ) cigcnvectois associated with the eigenvalues \i 2 = A + 1 A + 2 M by 
(noisi subspace cigcnvoctoi matiix) so that 

B = [B,B,] (2 42) 


Now fiom (2 42) and (2 39) since B^ B — I 
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[B B f{ARA^)[B Bn] + ct^I = Diag[Ai A^] (2 43) 

riio above equation leads to tlie following two equalities 

ARA^B = 5 Ai - S (o-a/) 

= B Diag[^i jii^] (2 44) 

ARA^B, = B,A2- Bnia2l)^0 (2 15) 

1 leis( note lliit m ( qns (2 44) and (2 45) The / s denote identity matrices of orders A and 
(M-A) Hspetlively Since AA is of lankA it follows from (2 15) that the (M x (ikf- A )) matrix 
A^^ B„ must satisfy 

A^^Bn = Q (2 46) 

Since BY B^ = I eqn (2 44) implies {Bf A)R{A^ Bs) = (Dta£()[^i /Ja Aa] whirh in tnin implies 
lhal the (A X A) matrix Bg is nonsmgulai Thus it follows from (2 44) that 

A = A,Diag[/q fJ.K]{A^ B )~^ R-^ 

= BgC (2 47) 

when C = ])aig[/Ji /ia]( IS nonsmgulai It ran be seen in the following section that 
tin columns o( matiix C are related to the eigenvectois of the matrix whose eigenvalues etimates 
the paiamc ters(in this case DOA’S) in the standard model of the signal correlation matrix 

2 4 2 Derivation of the Algorithm 
Aftei dehning, 

a(^) = [1 2 (2 4S) 
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it IS clesi Iioin (2 46) that /j — 2=12 K are the Mroots of the (M — l)th degire 

polynoiTinl equations 


[a{z)fb = bo + biz+ +bM-iz^~'^ = 0 (2 49) 

wheic b = [ 6 o Ih bM-if is any one of the columns of Theiefore 

A^b = 0 (2 50) 

Each column of i? ^ is associated with an equation similar to the eqn (2 50) The set of these 
2{AI — h) ccju itionstan be wiitten as 

[c/('')]‘[/;| &2(M-/v)] = 0 

and from (2 50) 


^2{M-A)] = 0 


(2 51) 


Ilcie it IS issumed that 6 m-i 7 ^ 0 in (2 49) for if 6 m-i = 0 the monomial of highest degree is 
selected Multiply (2 49) by 2 : dividing through out by 6 m~i and rearrange as follows 

z" + = 0 (2 52) 

bi\i-i Ojvf-i 

A(l( t issoci Uiug i companion maliix D of ordei M with the polynomial on the left hand side of 
(2 52) 


0 10 0 
0 0 1 0 


0 0 0 1 

Q ““6n —61 

^ bu-^l i'M-l ^M-1 


(2 53) 


it IS cle ir that the polynomial equation in (2 52) may be written as 
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det[z I — D] = Q The clesued roots z = 2=12 A on the unit circle are also the A 

unit maftirilude eigenvalues ^2 = 12 A of A> 

2 4 3 Theorem 1 

The (A X A ) matiiXjBjDB is nonsmgular whose eigen values uniquely estimate the A distinct 
paiameteis 2=12 /r in the standard model of the signal correlation niatiix 

Foi the proof please lefei [6] 



Chapter 3 

Evaluation Of C-R Bound For One 
Source Case 


These bounds apidy to unbiased cstimatois They are useful m that the best possible performance 
(smilhst vuniuc) loi on unbiased estimatoi may be determined [S] Estimators whose vaiiance is 
close oi ecpials the bound tan then said to be an optimal one 

3 1 Joint Likelyhood Function 

1 he inodtl foi diiection of ariival pioblem as we have shown in eep (2 3) is 

h 

m + (* - l)^r)) + « 

r=l 

Eoi OIK soutce case lei us define 


U) = OJi = Tf COS 


(3 2) 


Wf know fiom chap(ti2 of Vin tiecs [7] tint if 

x{t) = y(t,A) + n{f) (3 3) 

wheic n(L) is zero mean white Gaussian Noi&c with power spectral density ^,i:(f) is the observed 
signal y(l A) is the tiansmitted waveform and A is the unknown parameter winch has to be 
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estumtecl Then the loft likelihood Innction is given by 

A)cZi (34) 

T.0 IiTul 1 h( 111 1X1 in uni likoliliood estimate of A we have to hud the value of A for which ec{ii( 1 1 ) 
IS maximum A iiecessiiy but not sufficient condition can be obtained by differentiating eqn (3 4 ) 
With lespect to A and equating it to zero 

= 0 (3 5) 

Above c ((11 ilioii (iiovides ns tlu values of A al whirh the derivative of the log likelihood function is 

/( lO 


Now when wo have iiioie than one simultaneous observations as in our sensor array case we 
should find t Im loint log likelihood function We can ex;tend the eqn(3 4) directly to the case wheie 
oiu ui ly li IS N sdisois so that we have N obseivations Then 


dlnr^/gjXlA) 

dA 


2 

No 



k(t) -y{t A)] 


A) 

dA 


dt 


(3 6) 


wlicit (he unknown paiaineter and a.,(i) = y {t u) n (t) 0 < t < T where n (t) is 

/CIO mom while" C aussiaii noise with p s d ^ 


11 A IS 1 lion i uuknii puaindci (hen the ( lamei Rao bound requires that 


Var(f^-^ > -E[ 


d^lnP^u{X/A)-^ 


dA^ 


(3 7) 


Snbstilnting A is u) Now 




dHaP^UX/W)^ 2 




N i-T 




d\i-t W),^ [T dy^(tv^2. 




if (3S) 


wheic wo assume the clciivatives exist In the fust term 
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Ci'c.it) - y,{t w)] = E[n {t)] = 0 (3 9) 

In the second tenii tlioic aie no landoni quantities hence the expectation operation gives the 
integial itself 
Siibstilutnig in oqn (3 7) 


Var{u - w) > 


No 


1 01 tlio single sonice c xso 


(3 10) 


(/ (/ U) = VlPip{-{jWrt + (l - l)w)) 

winch implies - ‘ = -j(*- l)2/ieTp(— (juv^ + (* - l)w)) 

Substitnling in the eqn (3 10) 


01 


l '^ ai (cu — w) > 


No 




3JVo 

|i2TiV(iV- l)(2iV- 1) 


(311) 


(3 12) 


Var{u> ~ u))> 


6 

[TiV(iV - l)(2iV- l)sra] 


(3 13) 


«<'< ' 

2 

As expected, the ( 11 bound loi the frequency estimator is inversly proportional to snr (7^) 
rncl N (nnmbci of scnsois) 


Poi single souice case with M — 3 and M ~ 5 the perfoimance of MUSIC estimator has been 
sunnKtecl loi various sni ranges from 5 to 30dB 
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Chapter 4 


Simulation Results 

In MUSIC iiiothod it lequiies compulition of perfect oithogonal eigenvectors correspond to inin 
iinuin cigdiviliio wlmh his inullipliiity of M — A To compute the eigenvalues/vectois for a 
lldiniti 111 111 iliiJv il H c|uiKs pi icticilly infinite nuinbei of plane rotations and leflections In x>iin 
ti|)l( lU 1 ilioiis IK stojiped when the off diagonal elements are negligible to the woiking accuracy 
[^] 


I uobi itciitivc algoiithin ind Householder method produce almost exact orthogonal eigenvec 
(ois but m ly be in k (ui ito when the matrix undei consideration has close eigenvalues [9][11] Also 
Ilouselioldci mol hod is lolitively faster than the latter one 

Ilencc no xlgoulhm is best legarding the computation of eigen values/vectors of complex ma 
tiix[10] 

I h( dime iision of ( he sign il snbspad is estimated by comparing the eigenvalues of the correla 
tion iin(ii\ ind using i lluoshold It can be seen fiom the plots shown in figs (4 1) and (4 2) that 
the two smilloi dgeu viluos ((oiiesponding to noise subspace) are not same Also from fig (4 2) 
wo ( an say th it R ink ol the imtiix So = 5 - (T 2 / is no longer equal to A and it is gieater than A 
riiG plots wcie shown foi vaiious sni’s 10, 20 lOdB And the angle of arrivals were at 30 60deg 

As wo h ivc seen Diioct Noise Subspace Basis (DNSB) method is simple and computationally 
elhcionl if we cm estimate the sensor noise variance 
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Sun Illations weic cairied out foi one source and two sources with sensor sizes of M = 3 4 and 5 
In all the cases the data length was 100 and the input sinusoidal frequencies were 100 and 250 Hz 
loi two source case and 100 Hz was for one source case Random numbers were taken foi arbitrary 
seed value 

The DNSB method perfoiming well for M == 3 with A = 1 and A* = 2 Even for M = 4 and 5 
it helps to resolve the peal s at coirect locations but spread is more relatively with MUSIC method 
This is because foi M > A + 1 the condition on rank of 5“ < A is not achieved for real data 
enviiomnenl 

4 1 Single Target Case(A =1) 

MUSIC method foi direction of arrival computes the function 

PH = (4 1) 

2^ =A+1 

will coiiespond to the true diiection of anival DNSB method resolves the peak by solving the 
eci nation 


wlieio So = ( S — cT^J) which implies 


1 

5o vi 
0 


= 0 


(4 2) 


S [tJl] = - s 


(43) 


whoie S~ IS tlu second column of sq matiix Above eciuation is solved for vi by using left inverse 
(S-)+ =:(5 -"s-)“^H^ of S'- 

Tlie noise eigenvectors V[ V 2 using the light shift operator are 


1 


■ 0 

Hi 


1 

0 


. '"1 - 
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lespectively 

Now by using oqn (4 1) we compute the powei spectrum The plot for single binusoicl (arrival angle 
GOdcg ) IS shown in fig (13) foi M = 3 It can be seen that the spread is less in case of DNSB 
method Ilian MUSIC method The snr was 20dB foi the sinusoid 

4 2 Double Target Case(A = 2) 

For M = 3 


Ul 

Vl 

wluie S' IS the M X K imtiix consists of second and thud columns of matrix 


(4 4) 


Due to the 1 ank degeneiacy in S~ matiix we have done patiial pivoting before finding out 
Noise subspice basis vectois iq We have searched foi laigest element in magnitude in the 

hisl column then by inteichangmg the rows we have put that element at (1 1) position Similarly 
second (olumn is seaiched foi the laigest element in magnitude and bringing that element at the 
(oiiesponcling (2 2) position Hence so on and so fourth Also the diagonally decreasing ordci in 
magnitude is assured Finally we have forced the M — K rows to zero 

By using left mveise we compute vi V 2 then noise eigenvector Fi 


Similatly foi M = 4 and M = 5 we hnd out noise eigenvectors from that the eqn (4 1) to resolve 
I Ik peaks 


Foi M = lanclA = 2 MUSIC method fails to resolve the peaks at correct locations and the 
spiead IS moic also Hence in this case DNSB method is the ideal choice Further more foi 
M = 4 5anclA = 2 none of the Iwo methods able to resolve the peaks Only one of the two angles 
could able lo detect by them ihese plots are shown m figs (4 4), (4 6) and (4 6) 


Angle clistiibiition plots were shown in figs (4 7) to (4 12) The SNR was 20dB for each sinusoid 
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Clculy fiom these plots we can say that the DNSB is the best estimator for M = 3andA = 2 as 
it has lowei vuniicc we took 100 iteiations to plot these figs From the plots variance and bias 
were found using the statistical formnlae The two arrival angles were at 30 and 60 degs But for 
M = 5uiclA = 2 the same cannot conclude as the lank of the matrix 9o = (5 - a^I) is not equal 
to K 

The two inetliods (MUSIC and DNSB) compaied on the basis of comutational time In case 
of DNSB method the cpn time found out for the noise subspace basis vi vk through Moore 
Peniose inverse fiom the niatiix And in MUSIC method it was to find out eigen values and 
eigen vcctois fiom the c on elation niatiix 5 

1 oi two soul CO ( ISO with scnsoi si/es 5 7 and 10 epu tune found out for the two methods 
MUSK til os II Ims, 24ins and lOms, whoieas DNSB takes 6 66ms 6 Sms and 7 "ims respectively 

The CR hound as given by the eqn (3 13) is plotted assumming the estimator is unbiased so 
that the mein sqimio cnoi is equal to the variance The estimator nearly attains the CR bound 
ibove some thieshold SNR (aboni lOdB for M = and about 20dB for M = 3) Unfortunately 
below the thieshold the mean sqnre enoi inci eases lapidly, which renders the estimator not useful 

Finally fiom the plots shown in (igs (4 13) and ( 1 14) the MUSIC algorithm gives acceptable 
pel foi mane c foi sin s ibovc lOdB 



Chapter 5 

Conclusions And Suggestions For 
Further Work 


We have iiiade an a,1;1;empt to estimate the DOA’s using the Eigenstnictnre based methods MUSIC, 
Direct Noise Siibspace Basis and Subspace Based Approach. 

The signals are narrowband and therefore the amplitude and phase modulations are assumed 
constant as the wavefronts travel across the array. As signal travels from one sensor to another it 
undergoes a phase change. This phase change is a function of the angle of arrival. We have tried 
to estimate this angle of arrival with the above mentioned three methods in the presence of zero 
mean white gaussian noise with variance All the three methods were verified with real data. 

In chapter2 we have breifly described the three algorithms analytically. The DNSB method 
though it fails to perform that well for M > K + 1 ns it was for M = A' + 1 it is simple and 
computationally efficient . 


In chapters analytical solution for for one source case has been derived. The estimate' turned 
to be unbiased assuming high SNR. 

The possible directions in which further study in this area can be done are: 

For M > K + 1 one can find out the ways to make DNSB method to resolve the peaks with 


better resolution. 
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Analytical solution foi CR bound for two source case with more than three sensors can also be 
derived to verify with simulation results. 
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